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ABSTRACT 

Context. The ultra-compact dipping source XB 1916-053 has an orbital period of close to 50 min and a companion star with a very 
low mass (less than 0.1 M 0 ). The orbital period derivative of the source was estimated to be 1.5(3) x 10 -11 s/s through analysing the 
delays associated with the dip arrival times obtained from observations spanning 25 years, from 1978 to 2002. 

Aims. The known orbital period derivative is extremely large and can be explained by invoking an extreme, non-conservative mass 
transfer rate that is not easily justifiable. We extended the analysed data from 1978 to 2014, by spanning 37 years, to verify whether a 
larger sample of data can be fitted with a quadratic term or a different scenario has to be considered. 

Methods. We obtained 27 delays associated with the dip arrival times from data covering 37 years and used different models to fit the 
time delays with respect to a constant period model. 

Results. We find that the quadratic form alone does not fit the data. The data are well fitted using a sinusoidal term plus a quadratic 
function or, alternatively, with a series of sinusoidal terms that can be associated with a modulation of the dip arrival times due to 
the presence of a third body that has an elliptical orbit. We infer that for a conservative mass transfer scenario the modulation of the 
delays can be explained by invoking the presence of a third body with mass between 0.10-0.14 M 0 , orbital period around the X-ray 
binary system of close to 51 yr and an eccentricity of 0.28 ± 0.15. In a non-conservative mass transfer scenario we estimate that the 
fraction of matter yielded by the degenerate companion star and accreted onto the neutron star is /J = 0.08, the neutron star mass is 
> 2.2 M 0 , and the companion star mass is 0.028 M 0 . In this case, we explain the sinusoidal modulation of the delays by invoking the 
presence of a third body with orbital period of 26 yr and mass of 0.055 M 0 . 

Conclusions. From the analysis of the delays associated with the dip arrival times, we find that both in a conservative and non¬ 
conservative mass transfer scenario we have to invoke the presence of a third body to explain the observed sinusoidal modulation. We 
propose that XB 1916-053 forms a hierarchical triple system. 

Key words, stars: neutron - stars: individual (XB 1916-053) — X-rays: binaries — X-rays: stars - Astrometry and celestial me¬ 
chanics: ephemerides 


1. Introduction 


The X-ray source XB 1916-053 is a low-mass X-ray binary 
(LMXB) showing d ips and type-I X-ray bursts in its light curves. 
Using OSO 8 data, iBecker etafl (119771) observed type-I X-ray 
bursts, implying that the compact source in XB 1916-053 is a 
neutron star. Assuming that the peak luminosity of the X-ray 
bursts in XB 1916-053 is at the Eddington limit. [Stnale et alj 
(1988) derived a distance to the source of 8.4 kpc or 10.8 kpc, 
respectively, depending on whether the accreting matter has_cos- 
mic abundances or is extremely hydrogen-deficient. lYoshidal 
( 1993 ) inferred a distance to the source of 9.3 kpc studying 
the photospheric radius expansion of the X-ray bursts in XB 
1916-053 (see also iBarret et al]fl996h . XB 1916-053 was the 
first LMXB in which periodic absorptio n dips were detected 
(IWalter et al.ll 19821 IWhite & Swank! 19821) . These dips represent 
a decrease in the count rate in the light curve caused by periodic 
absorption of the X-ray emission produced in the inner region 
of the system. The photoelectric absorption occurs in a bulge at 
the outer radius of the accretion disc where the matter streaming 
from a companion star impacts. 


Accurate analysis of data sets from many X-ray satellites in 
the last 30 years have found different values for the X-ray pe¬ 
riod: IWalter et al.l (|198 2|) found a period close to 2 985 s, us¬ 
ing Einstein data; IWhite & Swankl (119821) e stimated a period o f 
3 003.6 + 1.8 s for the strongest dips, while lSmale et a l. (1989), 
analysing GINGA data, derived a period of 3 005.0 ± 6.6 s. 
IChurch et al.l (1997), analysing ASCA data, found an orbital pe¬ 
riod of 3 005 ± 10 s. The X-ray light curve of XB 1916-053 also 
shows secondary dips occurring approximately half a cycle away 
from the primary dips with a certain variability in phase (see 
lGrindlavil989l) . No eclipses were found; this constrains the or¬ 
bital inclination of the system between 60° and 80°. 


The optical counterpart of XB1916-053 was discovered by 
[Grindlav et al.l ( 19871). a star with a V magnitude of 21 already 
noted by Walter et all (1 19821). Usin g thermonuclear flash models 
of X-ray bursts. ISwank et all (119841) argued that the companion 
star is not hydrogen exhausted and suggested a_companion star 
mass of 0.1 M 0 . Furthermore. iPaczvnski & Sienkiewiczi d 19811) 
showed that X-ray binary systems with orbital periods shorter 
than 81 min cannot contain hydrogen-rich secondary stars. 
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Bin time: 64.00 s Bin time: 64.00 s 
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Fig. 1 . Chandra/LEG light curves of XB 1916-053 during the two observations performed in 2013, i.e. obsid. 15271 (left) and 15657 (right). The 
bin time is 64 s. A type-I X-ray burst that occurred during the obsid. 15271. 


A modulation in the optical light curve with a peri od of 
3 027.4 + 0.4 s was discovered bv iGrindlav et alj (Il988t) . The 
1% discrepancy between the optical and X-ray period of XB 
1916-053 was explained bv IGrindlav et al.l (1988) invoking the 
presence of a third body with a period of 2.5 d and a retrograde 
orbit that influences the matter streaming from the companion 
star. The same authors also suggested the alternative scenario 
in which the disc bulge precesses around the disc with a pro¬ 
grade period equivalent to the beat period between the optical 
and X-ray period. [Whitel ( 1989 ) suggested the possibility that 
a precessing elliptical disc exists in XB 1916-053, and that the 
variation in the projected area of this disc causes optical modula¬ 
tion. iGfikmaneCal] (1 1 9951) showed the st ability of the optical pe¬ 
riod over seven years. Chou et ah (2001), analysing Rossi X-ray 
Timing Explorer (RXTE) data taken in 1996, found several peri¬ 
odicities including one at 3 026.23 + 3.23 s, which was similar to 
the optical modulation at 3 027 s. The centroid of these peaks in 
the periodogram associated with the 3 000 s period implies that 
there is a modulation with a fundamental period close to 3.9 d, 
as already noticed bv iGrindlavl (119921) also in the optical band. 
The period of 3.9 d is interpreted as the beat period between the 
optical and X-ray periods. Furthermore. IChou et al.l ( 2001 ). fold¬ 
ing the RXTE light curves at the 3.9 d period, found changes in 
the dip shape following this modulation. Those authors also in¬ 
dicated that the dip-phase change, with a sinusoidal period of 
6.5 ± 1.1 d from Ginga 1 990 September observations (Yoshida 
il993l lYoshida et al.lll995l) . may be associa ted with the subhar¬ 
monic of the 3.9 d period. IRetter et all (120021) detected a further 
independent X-ray period at 2 979 s in the RXTE light curves 
of XB 1916-053, which was mistakenly identified by Chou et al. 
(2001) with a 3.9 d sideband of the 3 000 s period. IRetter et ah 
(2002) suggested that the period at 2 979 s could be explained 
as a negative super-hump assuming the 3 000 s period is the or¬ 
bital period with a corresponding beat period of 4.8 d. The same 
authors suggested that the 3.9 and 4.8 d periods could be the 
apsidal and nodal precession of the accretion disc, respectively. 

Finally, the source also showed along-term 198.6+1.72d pe¬ 
riodicity in X-rays (iPriedhorskv & Terre 111! 19841) . which has not 


been confirmed by further observations (see lRetter et al.l1200 2). 
To date the spin period of the neutron star in XB 1916-053 is not 
known. iGallowav et al.l (1200 ll) . analysing a Type-I X-ray burst, 
discovered a highly coherent oscillation drifting from 269.4 Hz 
up to 272 Hz. Interpreting the asymptotic frequency of the oscil¬ 
lation in terms of a decoupled surface burning layer, the neutron 
star could have a spin period around 3.7 ms. 

iHu et al.l (i2008l) inferred that P,„u !P„ r h = (1.62±0.34)x 10~ 7 
yr 1 by analysing archival X-ray data from 1978 to 2002 and 
adopting a quadratic ephemeris to fit the dip arrival times. In this 
work, we update the previously determined ephemeris using data 
from 1978 to 2014. We show that the quadratic ephemeris does 
not fit the dip arrival times and find that a sinusoidal component 
is necessary to fit the delays. We suggest the presence of a third 
body that influences the orbit of the X-ray binary system XB 
1916-053. 


2. Observations and data reduction 

We used all the available X-ray archival data of XB 1916-053 
to study the long-term change of its orbital period. The last 
ephemeris of the source was reported by iHu et al.l ( 2008 ) who 
used archival data from 1978 to 2002. We analysed more than 
37 years of observational data from 1978 to 2014. The data have 
been obtained from the HEASARC (NASA’s High Energy As¬ 
trophysics Science Archive Research Center) website and have 
been reduced using the standard procedures. In particular, we re¬ 
analysed the data used bv lHu et aD (120081 ). collected from 1998 
to 2002, and added new data spanning up to 2014 (see Tab. Q}. 
We obtained 27 points from all the analysed observations. The 
data collected by RXTE, Ginga, EXOSAT, Einstein, and OSO- 
8 were downloaded from HEASARC in light-curve format. We 
used the standard-1 RXTE/PCA background-subtracted light 
curves, which include all the energy channels and have a time 
resolution of 0.125 s. All the pointing observations were used 
except for P70034-02-01-01, P70034-02-01 -00, and P93447-01 - 
01-00 due to the absence of dips in the corresponding light 
curves. The EXOSAT/ME light curves cover the energy range 
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Table 1. Observation Log. 


Point 

Satellite/Instrument 

Observation 

Start Time 
(UT) 

Stop Time 
(UT) 

Tfold 

(MJD,TDB) 

1 

OSO-8/GCXSE 


1978 Apr 07 21:16:05 

1978 Apr 14 22:20:37 

43 609.408575724435 

2 

Einstein/IPC 


1979 Oct 22 04:52:01 

1979 Oct 22 06:58:30 

44 168.24670380917 

3 

Einstein/IPC 


1980 Oct 11 04:08:51 

1980 Oct 11 09:07:19 

44523.27644368849 

4 

EXOSAT/ME 


1983 Sep 17 15:07:25 

1983 Sep 17 21:29:49 

45 594.765324269885 

5 

EXOSAT/ME 


1985 May 24 12:26:21 

1985 May 24 21:30:23 

46 209.612747685185 

6 

EXOSAT/ME 


1985 Oct 13 13:53:16 

1985 Oct 13 22:34:04 

46351.75944524423 

7 

Ginga/LAC 


1988 Sep 09 15:47:56 

1988 Sep 10 16:01:16 

47414.165911835925 

8 

Ginga/LAC 


1990 Sep 11 15:04:35 

1990 Sep 13 09:18:11 

48 146.51075733274 

9 

ROSAT/PSPC 

RP400274N00 

1992 Oct 17 13:05:47 

1992 Oct 19 15:24:20 

48 913.59379352164 

10 

ASCA/GIS3 

40004000 

1993 May 02 18:11:00 

1993 May 03 09:46:17 

49110.082393510115 

11 

RXTE/PCA 

P10109-01 -01-00, P10109-01-02-00, 

P10109-01 -04-01, P10109-01-04-00, 

P10109-02-01 -00, P10109-02-02-00, 

P10109-02-03-00, P10109-02-04-00, 

P10109-02-05-00, P10109-02-06-00, 

P10109-02-07-00,P10109-02-08-00, 

P10109-02-09-00, P10109-02-10-00, 
P10109-02-10-02 

1996 Feb 02 00:14:56 

1996 May 23 11:20:00 

50174.74129123185 

12 

RXTE/PCA 

P10109-01 -05-00, P10109-01 -06-00, 

P10109-01 -07-00, P10109-01 -08-00, 
P10109-01-09-00 

1996 Jun 01 17:38:40 

1996 Oct 29 11:00:34 

50310.596956288645 

13 

BeppoSAX/MECS 

20106001 

1997 Apr 27 21:00:06 

1997 Apr 28 19:51:02 

50566.35264963594 

14 

RXTE/PCA 

P30066-01 -01 -04, P30066-01 -01 -00, 
P30066-01-01-01, P30066-01 -01 -02, 
P30066-01 -01 -03, P30066-01 -02-00, 
P30066-01 -02-01, P30066-01 -02-02, 
P30066-01-02-03 

1998 Jun 23 23:06:40 

1998 Jul 20 15:35:55 

51001.306447481845 

15 

RXTE/PCA 

P30066-01 -02-04, P30066-01 -02-07, 
P30066-01 -02-08, P30066-01 -03-00, 
P30066-01 -03-01, P30066-01 -03-02, 
P30066-01 -03-03, P30066-01 -03-04, 
P30066-01-03-05, P30066-01-04-00 

1998 Jul 21 07:11:44 

1998 Sep 16 02:52:32 

51043.70980975036 

16 

RXTE/PCA 

P30066-01 -05-01, P30066-01 -05-00, 
P30066-01 -06-00, P30066-01 -06-01, 
P30066-01-07-00, P30066-01-07-01 

2001 May 27 08:14:47 

2001 Jul 01 19:15:33 

52074.07302734295 

17 

BeppoSAX/MECS 

21373002 

2001 Oct 01 03:40:16 

2001 Oct 02 07:01:06 

52 183.72270184033 

18 

RXTE/PCA 

P50026-03-01-00, P50026-03-01-01 

2001 Oct 01 10:35:44 

2001 Oct 01 22:16:03 

52 183.684644754605 

19 

RXTE/PCA 

P70034-02-02-01, P70034-02-02-00 

2002 Sep 25 00:43:12 

2002 Sep 25 09:31:12 

52542.21332826887 

20 

XMM/Epic-pn 

0085290301 

2002 Sep 25 04:18:29 

2002 Sep 25 08:28:27 

52542.266295747205 

21 

INTEGRAL/JEM-X 


2003 Nov 09 09:04:11 

2003 Nov 20 12:18:01 

52957.945226848465 

22 

Chandra/HETGS 

4584 

2004 Aug 07 02:34:45 

2004 Aug 07 16:14:53 

53 224.59478392645 

23 

Suzaku/XISO 

401095010 

2006 Nov 08 06:09:51 

2006 Nov 09 02:42:02 

54048.3655207864 

24 

RXTE/PCA 

P95093-01-01-00, P95093-01-01-01 

2010 Jun 19 13:41:52 

2010 Jun 21 07:21:46 

55 367.43875650959 

25 

Chandra/LETGS 

15271, 15657 

2013 Jun 15 13:56:17 

2013 Jun 18 05:13:17 

56459.89915961875 

26 

Swift/XRT 

00033336001 

2014 Jul 15 08:04:57 

2014 Jul 15 22:36:46 

56853.63959388178 

27 

Suzaku/XISO 

409032010, 409032020 

2014 Oct 14 16:49:56 

2014 Oct 22 2:40:56 

56949.56345974802 


between 1 and 8 keV and have a time bin of 16 s. The Ginga/LAC 
light curves cover the 2-17 keV energy band. We only used the 
data from the top layer and the light curves binned at 16 s. We 
downloaded the ROSAT/PSPC events, and extracted the corre¬ 
sponding light curve using the FTOOLS xselect. The Medium 
Energy Concentrator Spectrometer (MECS) onboard the Bep- 
poSAX satellite observed XB 1916-053 two times, in 1997 Apr 
27-28 and 2001 Oct 01-02. Using xselect, we extracted the 
source light curves from a circular region centred on the source 
and with a radius of 4', no energy filter was applied to the data. 
The BeppoSAX/MECS light curves were obtained using a bin 
time of 2 sec. ASCA observed XB 1916-053 in 1993 May 02- 
03; we used the events collected by the GIS3 working in medium 
bit rate to extract the corresponding light curve. The OSO-8 light 
curve was obtained using the combined observation of the B and 
C detectors of the GSFC Cosmic X-ray Spectroscopy experiment 
(GCXSE). The light curve covers the 2-60 keV energy range. 
The Einstein light curve was obtained from events collected by 
the Image Proportional Counter (IPC) in the 0.2-3.5 keV energy 
range. 

We applied barycentre corrections to the whole data set 
adopting the source position of XB 1916-053 shown by 
llaria et all d2006i) . For the RXTE/PCA light curves we used the 
ftools faxbary. The barycentre corrections for the ASCA and 


ROSAT data were obtained using the ftool timeconv and the 
tool bct+abc, respectively. All the other data sets were corrected 
using the ftool earth2sun. Finally, we excluded the time inter¬ 
vals containing X-ray bursts from each analysed light curve. 

The Chandra satellite observed XB 1916-053 three times. 
The first time was on 2004 Aug 07 from 2:34:45 to 16:14:53 UT 
(obsid 4584). The observation had a total integration time of 50 
ks and was performed in timed graded mode. The spectroscopic 
analysis of this data set was discussed bv llaria et al.l (1200 6). We 
reprocessed the data and applied the barycentre corrections to the 
event-2 file using the Chandra Interactive Analysis of Observa¬ 
tions (CIAO) tool axbary. In addition, we extracted the summed 
first-order medium energy grating (MEG) and high energy grat¬ 
ing (HEG) light curves filtered in the 0.5-10 keV energy range 
using the CIAO tool dmextract. The last two Chandra observa¬ 
tions of XB 1916-053 (obsid 15271 and 15657) were performed 
between 2013 June 15 13:56 and June 18 5:13 UT and have ex¬ 
posure times of 60 and 30 ks, respectively. We reprocessed the 
data and applied the barycentre corrections to the event-2 file us¬ 
ing axbary. Moreover, we extracted the first-order low energy 
grating (LEG) light curve in the 0.5-5 keV energy range using 
dmextract. We show the Chandra/LEG light curve in Fig. Q] 
Very intense dipping activity is present during the two observa¬ 
tions. A type-I burst occurred during the obsid. 15271. 
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Table 2. Best-fit parameters obtained fitting the dips in the folded light curves. 


Point 

Phase Interval 

C, 

count s 1 

c 2 

count s -1 

C 3 

count s 1 

01 

02 

03 

04 

xi d (.d.o.f.) 

1 

0.7-1.7 

7.23 ± 0.06 

6.01 ±0.09 

7.59 ± 0.06 

1.086 ±0.012 

1.113 ±0.012 

1.290 ±0.016 

1.388 ±0.015 

1.64(193) 

2 

0.7-1.7 

5.74 ± 0.06 

2.94 ±0.13 

5.75 ±0.11 

1.133 ±0.004 

1.170 ±0.004 

1.217 ±0.005 

1.263 ±0.006 

1.44(152) 

3 

0.1-1 

11.24 ±0.08 

8.48 ±0.12 

11.14 ±0.09 

0.416 ±0.006 

0.450 ± 0.005 

0.574 ± 0.006 

0.603 ± 0.008 

1.90(194) 

4 

0.8-1.8 

23.02 ±0.12 

13.7 ±0.4 

23.73 ±0.11 

1.203 ±0.003 

1.239 ±0.004 

1.277 ±0.003 

1.321 ±0.003 

1.91(194) 

5 

0.095-0.8 

14.6 ±0.4 

8.9 ±0.3 

17.7 ±0.3 

0.251 ±0.015 

0.329 ±0.015 

0.500 ± 0.008 

0.576 ± 0.008 

4.70(114) 

6 

0.1-0.8 

25.6 ±0.2 

20.3 ± 0.4 

27.1 ±0.2 

0.454 ± 0.007 

0.487 ± 0.008 

0.561 ±0.007 

0.605 ± 0.007 

3.03(133) 

7 

0.5-1.1 

72.9 ± 1.1 

42.8 ± 1.0 

75.6 ± 1.0 

0.658 ± 0.004 

0.680 ± 0.004 

0.840 ± 0.009 

0.940 ±0.010 

16.9(84) 

8 

0.5-1.2 

107.1 ± 1.0 

60.6 ± 1.3 

106.3 ±0.8 

0.657 ± 0.007 

0.805 ± 0.008 

0.903 ± 0.004 

0.954 ± 0.005 

47.8(138) 

9 

0-1 

6.05 ± 0.07 

1.60 ±0.06 

6.15 ±0.07 

0.340 ± 0.008 

0.535 ± 0.006 

0.640 ± 0.004 

0.752 ± 0.005 

4.77(294) 

10 

0.6-1.4 

9.3 ±0.3 

0.85 ± 0.05 

9.3 ±0.3 

0.765 ± 0.004 

0.809 ± 0.002 

1.022 ±0.003 

1.100 ±0.006 

12.08(234) 

11 

0-1 

41.00 ±0.11 

31.35 ±0.11 

45.06 ± 0.08 

0.207 ± 0.004 

0.323 ± 0.004 

0.499 ± 0.002 

0.611 ±0.003 

12.61(506) 

12 

0.1-1.1 

36.60 ±0.10 

21.69 ±0.13 

37.45 ±0.11 

0.443 ± 0.003 

0.569 ± 0.003 

0.710 ±0.002 

0.808 ± 0.002 

11.33(505) 

13 

0.95-1.95 

0.986 ±0.013 

0.035 ± 0.003 

0.982 ±0.014 

1.238 ±0.004 

1.340 ±0.002 

1.537 ±0.002 

1.647 ±0.005 

2.14(249) 

14 

0.24-0.75 

27.6 ± 0.2 

14.55 ±0.10 

27.7 ± 0.2 

0.313 ±0.003 

0.433 ± 0.002 

0.582 ± 0.002 

0.705 ± 0.003 

11.10(255) 

15 

0.15-1 

36.51 ±0.10 

25.3 ±0.2 

38.31 ±0.08 

0.381 ±0.004 

0.566 ± 0.005 

0.598 ± 0.003 

0.738 ± 0.003 

10.43(420) 

16 

0.35-0.9 

25.4 ± 0.2 

17.52 ±0.10 

26.9 ± 0.2 

0.419 ±0.004 

0.489 ± 0.004 

0.697 ± 0.003 

0.763 ± 0.003 

1.90(274) 

17 

0.9-1.9 

1.030 ±0.009 

0.31 ±0.02 

1.016 ±0.007 

1.045 ±0.004 

1.160 ±0.004 

1.178 ±0.004 

1.272 ±0.004 

1.11(249) 

18 

0-1 

23.75 ±0.14 

11.4 ±0.3 

24.40 ±0.13 

1.092 ±0.002 

1.113 ±0.002 

1.168 ±0.005 

1.281 ±0.005 

1.57(505) 

19 

0.6-1.6 

28.8 ± 0.2 

20.1 ±0.2 

29.8 ± 0.2 

0.953 ± 0.004 

0.985 ± 0.004 

1.171 ±0.006 

1.249 ±0.006 

3.06(505) 

20 

0.1-1.1 

69.9 ± 0.5 

30.4 ± 0.5 

72.4 ± 0.5 

0.406 ± 0.004 

0.506 ± 0.003 

0.671 ±0.003 

0.730 ± 0.003 

17.84(505) 

21 

0.2-1.2 

0.0284 ± 0.0002 

0.0256 ± 0.0004 

0.0284 ± 0.0002 

0.57 ± 0.02 

0.63 ± 0.02 

0.708 ± 0.014 

0.733 ±0.011 

0.766(144) 

22 

0.3-1.3 

9.92 ± 0.03 

8.5 ±0.2 

10.04 ± 0.03 

0.796 ± 0.005 

0.835 ± 0.007 

0.844 ±0.013 

0.923 ± 0.007 

1.71(144) 

23 

0.9-1.9 

15.10 ±0.05 

11.4 (fixed) 

15.61 ±0.05 

1.263 ±0.006 

1.38 (fixed) 

1.38 (fixed) 

1.521 ±0.006 

2.41(505) 

24 

0.9-1.9 

30.18 ±0.13 

23.1 ±0.3 

20.49 ±0.15 

1.322 ±0.003 

1.344 ±0.003 

1.499 ±0.002 

1.451 ±0.002 

1.22(505) 

25 

0.09-0.65 

2.52 ± 0.03 

0.759 ±0.010 

2.69 ± 0.04 

0.194 ±0.003 

0.314 ±0.002 

0.525 ± 0.002 

0.598 ± 0.003 

2.39(567) 

26 

0.6-1.6 

12.61 ±0.11 

4.4 ± 0.6 

12.5 ±0.2 

1.025 ±0.010 

1.27 ±0.02 

1.295 ±0.004 

1.307 ±0.004 

4.08(171) 

27 

0.84-1.6 

5.08 ± 0.04 

1.607 ±0.015 

5.18 ±0.02 

0.909 ± 0.003 

1.089 ±0.002 

1.2080 ±0.0011 

1.2904 ±0.0014 

3.29(382) 


Bin time: 16 s 



Fig. 2. Suzaku/XISO light curve of XB 1916-053 during the long obser¬ 
vation on 2014 Oct. 


The X-ray Multi-Mirror Mission-Newton (XMM-Newton) 
observed XB 1916-053 on 2002 Sep 25 from 3:55 to 8:31 UT 
and the European Photon Imaging Camera (EPIC-pn) collected 
data, in timing mode, over ~ 17 ks of exposure. An extensive 
study of this observation was performed bv lBoirin et al.l (120041) . 
We reprocessed the data, extracted the 0.5-10 keV light curve, 
and applied barycentre corrections to the times of the EPIC-pn 
events with the Science Analysis Software (SAS) tool barycen. 

Suzaku observed XB 1916-053 twice, the first time on 2006 
Nov 8 (obsid. 401095010) and the second time from 2014 Oct 14 
to 22 (obsid. 409032010 and 409032020). The first observation 


Table 3. Journal of the X-ray dip arrival times of XB 1916-053. 


Point 

Dip Time 
(MJD;TDB) 

Cycle 

Delay 

(s) 

1 

43 609.4168(12) 

-187 551 

772 ± 74 

2 

44 168.2535(5) 

-171460 

792 ± 28 

3 

44523.2941(5) 

-161237 

641+42 

4 

45 594.7744(3) 

-130385 

449+ 18 

5 

46 209.6271(13) 

-112681 

193 + 112 

6 

46 351.7778(9) 

-108 588 

352 ± 52 

7 

47 414.193(2) 

-77 997 

162+132 

8 

48 146.539(3) 

-56 910 

47+182 

9 

48 913.6127(10) 

-34 823 

-140 + 59 

10 

49109.1148(12) 

-29 165 

-48 ± 76 

11 

50 174.7555(5) 

1490 

-50 ± 46 

12 

50 310.6187(4) 

5 402 

-17 + 37 

13 

50566.3680(4) 

12 766 

-69 ± 39 

14 

51001.3241(5) 

25 290 

-15+40 

15 

51043.7292(5) 

26511 

-9+45 

16 

52 074.0935(3) 

56179 

151+29 

17 

52 183.7349(3) 

59 336 

107 ± 28 

18 

52 183.7008(2) 

59 335 

162+19 

19 

52 542.2168(4) 

69 658 

227 ± 39 

20 

52 542.2860(11) 

69 660 

202 ± 98 

21 

52 957.9679(8) 

81629 

327 ± 69 

22 

53 224.6246(4) 

89 307 

467 ± 34 

23 

54048.3791(5) 

113 026 

411+39 

24 

55 367.45218(15) 

151007 

593 + 13 

25 

56 459.9129(3) 

182463 

721+20 

26 

56 853.6454(8) 

193 800 

821+67 

27 

56 949.84670(10) 

196 570 

814 + 8 


Note — Epoch of reference 50123.00873 MJD, orbital period 
3 000.6511 s. 

has already been analysed by [Zhang et al.l d2014l) . while a study 
of the second observation has not been published yet. For both 
observations, we extracted the X-ray Imaging Spectrometer 0 
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(XISO) events from a circular region centred on the source and 
with a radius of 130". We applied the barycentre corrections to 
the events with the Suzaku tool aebarycen. We do not show the 
light curve of the first Suzaku observation since it was already 
shown bv lZhang et al.ld2014|) (Fig. 1 in their paper), however, we 
show in Fig.[2]the XISO light curve of the observation performed 
in 2014 Oct. The light curve indicates that a bursting activity is 
present in the first 200 ks of the observation and the persistent 
count rate decreases from 20 to 10 counts s' 1 . In the second part 
of the observation, the persistent count rate is quite constant at 7 
counts s _1 and an intense dipping activity is present. For the aim 
of this work, we selected and used the events from 250 ks to the 
end of the observation. 

Swift/XRT data were obtained as target of opportunity obser¬ 
vations performed on 2014 Jul 15 from 07:55:53 to 22:27:58 UT 
(ObsID 00033336001) for a total on-source exposure of ~ 6.3 
ks and on 2014 Jul 21 from 07:32:00 to 16:11:5 UT (ObsID 
00033336002) for a total on-source exposure of ~ 9.0 ks. The 
count rate in the first observation reaches 15 counts s -1 , with 
a mean at about 10 counts s -1 , due to the dips seen down to 
2 counts s 1 ; the second observation shows no dips and has a 
mean count rate of 7 counts s 1 . Since the data from ObsID 
00033336002 do not show dips we only used the first observa¬ 
tion in our analysis. The XRT data were processed with standard 
procedures (xrtpipeline v0.13.1), and with standard filtering 
and screening criteria with FTOOLS (v6.16). Source events (se¬ 
lected in grades 0 - 2 ) were accumulated within a circular region 
with a radius of 20 pixels (1 pixel ~ 2.36"). For our timing anal¬ 
ysis, we also converted the event arrival times to the solar system 
barycentre with barycorr. 

We selected a public data set of INTErnational Gamma-Ray 
Astrophysics Laboratory (INTEGRAL IWinkler et alJl2003h ob¬ 
servations performed in staring mode on XB 1916-053. Then, 
we analysed the data collected by the X-ray telescope JEM- 
X2 (iLund et al.ll2003l) . A total amount of 87 pointings (the to¬ 
tal observation elapsed time is ~ 310 ks) covered the INTE¬ 
GRAL revolutions 131, 133, and 134, which were carried out on 
2003 November 9-20. We performed the JEM-X2 data analysis 
using standard procedures within the Offline Science Analysis 
softw are (OSA10.0) distributed by the ISDC dCourvoisier et alJ 
120031) . We extracted the light curves with a 16 seconds bin- 
size in the energy range 3-10 keV, and after that we applied the 
barycentre corrections to the events using the tool barycent. 


3. Data analysis 

We analysed 27 light curves and folded the barycentric-corrected 
light curves using a trial time of reference and orbital period, 
Tf 0 id and Pq, respectively. For each light curve, the value of Tf„i d 
is defined as the average value between the corresponding start 
and stop time. We fitted the dips with a simple model consist¬ 
ing of a step-and-ramp function, where the count rates before, 
during, and after the dip are constant and the intensity changes 
linearly during the dip transitions. This model involves seven pa¬ 
rameters: the count rate before, during, and after the dip, called 
Ci, C 2 , and C 3 , respectively; the phases of the start and stop time 
of the ingress ( <p\ and fc), and, finally, the phases of the start and 
stop time of the egress (</>3 and (/> 4 ). The phase corresponding to 
the dip arrival time tpdip is estimated as (f>di P = (<p4 + 0i)/2. The 
corresponding dip arrival time is given by tdi p = T f a id + (f>di P Po- 
To be more conservative, we scaled the error associated with (f>di P 

by the factor Jx 2 ed to take a value of x 2 red °f the best-fit model 
larger than 1 into account. To obtain the delays with respect to 


a constant period reference, we used the values of the period 
Po = 3 000.651 l_s_and reference epoch Tq = 50 123.00873 MJD 
reported in lHu et akl (20081 ). We show the values of Tf„ui in Tab. 
Q] The best-fit parameters of the step-and-ramp function and the 
corresponding x 2 red are shown in Tab. [2] The inferred delays, in 
units of seconds, of the dip arrival times with respect to a con¬ 
stant orbital period are reported in Tab. [3] For each point we 
computed the corresponding cycle and the dip arrival time in 
days with respect to the adopted To. We show the delays vs. time 
in Fig.[3](left panel). 

Initially we fitted the delays with a quadratic function 
y(t) — a + bt + ct 2 , 

where t is the time in days with respect to 7o, a = ATq is the 
correction to To in units of seconds, b = AP/Po in units of s 
d 1 with A P the correction to the orbital period, and finally, c = 
1/2 P/Po in units of s d 2 , with P, that is the orbital period 
derivative. The quadratic form does not fit the data, we obtained 
^ 2 (d.o.f.) of 194.6(24). Here, and in the following, we scaled the 

uncertainties in the parameters by a factor Jx 2 e j to take a value 

of *2 of the best-fit model larger than 1 into account. The best- 
fit parameters are shown in the second column of Tab. [4] The 
corresponding quadratic ephemeris (hereafter LQ ephemeris) is 

3 000.65094(14) 

T dip (N) = MJD(TDB) 50 123.0096(3) +-- -N+ 

+2.37(12) x 10^ 13 A 2 , 

( 1 ) 

where N is the number of cycles, 50,123.0096(3) MJD is the 
new Epoch of reference, the revised orbital period is P — 
3,000.65094(14) s, and the orbital period derivative obtained 
from the quadratic term is P = 1.36(7) x 10 11 s/s. The obtained 
quadratic ephemeris is compatible with that reported bv lHu et all 
(2008). We show the best-fit curve in Fig. [3] (left panel) and the 
corresponding residuals in units of seconds in Fig. |3](right panel, 
upper plot). 

As we obtained a large value of the ^ 2 , we fitted the delays 
vs. time adding a cubic term to the previous parabolic function, 
i.e. 

y(t ) = a + bt + ct 2 + dt 2 , 

where a, b and c are above defined whilst the cubic term, d, is 
defined as P/(6Pq), and P indicates the temporal derivative of 
the orbital period derivative. Fitting with a cubic function, we 
obtained a^ 2 (d.o.f.) of 92.4(23) with a Ax 2 of 101.2 and an F- 
test probability of chance improvement of 4.2 xl0~ 5 with respect 
to the quadratic form. The best-fit values are shown in the third 
column of Tab. |4] The corresponding ephemeris (hereafter LQC 
ephemeris) is 


T dip (N) = MJD(TDB) 50 123.00870!“;™“+ 

3 000.65239(3) u „ , 

+-——————-A + 2.97(12) x 10- 13 A 2 - 2.2(4) x 10~ 22 A 3 ; 

86 400 

( 2 ) 

in this case we find an orbital period derivative of 1.71(7) x 10~ n 
s/s and its derivative is P = -3.8(0.7) x 10~ 2 ° s/s 2 . 

We also fitted the delays using a linear plus a sinusoidal func¬ 
tion having the following terms 


y(t) — a + bt + A sin 


2n 


Pmod 


(t - tf) 


(3) 
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Fig. 3. Left panel: dips’s arrival time delays versus time. The magenta, blue, black , and green curves are the best-fit curves obtained using the 
linear+quadratic (LQ), linear+sinusoidal (LS), linear+quadratic+sinusoidal (LQS), and linear+sinusoidal function taking into account a possible 
eccentricity (LSe), respectively. Right panel: observed minus calculated delays in units of seconds. The residuals, from the top to the bottom, 
correspond to the LQC, LS, LQS, and LSe function, respectively. 


Table 4. Best-fit values of the parameters of the functions used to fit the delays. 


Parameters 

LQ 

LQC 

LS 

LQS 


LSe 


a (s) 

78 + 23 

-2 1 +2A 

• -11.2 

584+ 157 

16 + 22 

229 ± 336 

56 ± 322 

1.1 ±299.2 

b(x 10~ 3 sd~‘ ) 

-4 + 4 

37.1 ±0.4 

-43 ± 23 

-4 + 3 

3+20 

3 + 19 

5+22 

c (xl(L 5 s d~ 2 ) 

1.70 + 0.09 

2.13 + 0.03 

- 

1.79 + 0.09 

- 

- 

- 

d (xl(L 9 s d~ 3 ) 

- 

-1.35+0.12 

- 

- 

- 

- 

- 

A(s) 

- 

- 

658 + 206 

130+15 

519 + 47 

548 ± 43 

577 ± 43 

U (d) 

- 

- 

3 897 + 332 

1 356 ± 203 

-3 723 + 1 100 

-3 150+1 116 

-2 923 ± 1 034 

Pmod (d) 

- 

- 

20409 + 3 381 

9 302 ± 752 

17 100 (fixed) 

18 600 (fixed) 

20 100 (fixed) 

m (deg) 

- 

- 

- 

- 

195 ± 26 

210 + 28 

217 + 27 

e 

- 

- 

- 

- 

0.26 ± 0.20 

0.28+0.15 

0.32 + 0.13 

T 2 (d-o.f.) 

194.6(24) 

92.4(23) 

63.7(22) 

39.4(21) 

51.8(21) 

48.2(21) 

45.8(21) 


Note — The reported errors are at 68% confidence level. The fit parameters of the delays are obtained using LQ (column 2), LQC (column 3), LS 
(column 4), LQS function (column 5), and LSe (columns 6, 7, and 8), respectively. 


where a and b are defined as above, A is the amplitude of the 
sinusoidal function in seconds, P mo d is the period of the sine 
function in days, and, finally, is the time in days referred to 
To at which the sinusoidal function is null. We obtained a value 
of^ 2 (d.o.f.) of 63.7(22) with a A ^ 2 of 131 with respect to the 
quadratic form. The best-fit parameters are shown in the fourth 
column of Tab. [4] The best-fit function is indicated with a blue 
curve in Fig. [3] (left panel) and the corresponding residuals are 
shown in Fig. [2 (right panel, the second plot from the top). The 
residuals are flatter than those obtained in the previous case. Us¬ 
ing the sinusoidal function, the dip time obtained from the OSO- 
8 observation is distant ~200 s from the expected value. The 
corresponding ephemeris (hereafter LS ephemeris) is 


T dip (N) = MJD(TDB) 50 123.01549(18)+ 


3 000.6496(8) 
86400 


IV + A sin 


2ji 

Nmod 


-N-< 


(4) 


where N mo d = Pmod/Po = 587 659.53 ± 97 351.67 and 0 = 
Int^/Pmod = with P, no d = 55.9 ± 9.3 yr. This functional form 
significantly improves the fit, even though it does not take the 
possible presence of an orbital period derivative into account. 

We added a quadratic term to take the possible presence of an 
orbital period derivative and fitted the delays into account, using 
the function 


y(t) = a + bt + cr + A sin 


2 7T 


Pmod 


U- 


(5) 


We obtained a value of^ 2 (d.o.f.) of 39.4(21) and a F-test prob¬ 
ability of chance improvement with respect to the LS ephemeris 
of 1.7 x 10~ 3 . The best-fit parameters are shown in the fifth 
column of Tab. [4] The best-fit function is indicated with a 
black curve in Fig. [3] (left panel) and the corresponding resid¬ 
uals are shown in Fig. [3] (right panel, the third plot from the 
top). The corresponding linear+quadratic+sinusoidal ephemeris 


Article number, page 6 of[[3] 




















R. Iaria et al.: Is XB 1916-053 part of hierarchical triple system? 



Fig. 4. Folded RXTE/ASM light curve of XB 1916-053 in the 3-5 and 5-12.2 keV energy range (top and middle panels). The corresponding 
hardness ratios_(HRs) are plotted in the bottom panels. The left and right plots show the folded RXTE/ASM light curve using the ephemeris 
discussed bv lHu et aL ( 2008 ) and LQ ephemeris (eq.|T} shown in the Sect. 3, respectively. Each phase-bin is about 50 s. 


Folded period: 3000.65239 s 
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Fig. 5. Left and right plots show the folded RXTE/ASM light curve using LQC ephemeris (eq. □ and LS ephemeris (eq. 0. respectively. Each 
phase-bin is about 50 s. 




0 0.5 1 1.5 2 


Orbital Phase 


(hereafter LQS ephemeris) is 


Tdi P (N) = MJD(TDB) 50 123.0089(3) + 


3 000.65126(10) 


+2.50(12) x 10~ 13 /V 2 + A sin 


86400 
2n 


N+ 


N m od 


N-<t> 


( 6 ) 


with N mo d = 267 837.87 ± 21 652.90 and <f> = 0.92 ± 0.16. The 
corresponding orbital period derivative is P = 1.44(7) x 10 11 s/s 
and the period of the modulation is P mo d = 25.5 + 2.1 yr. 

Our analysis of the delays suggests that a quadratic or a 
quadratic plus a cubic term do not fit the delays. A better fit is ob¬ 


tained using a sinusoidal function with a period close to 20 000 
d and, finally, adopting a sinusoidal plus a quadratic term, we 
obtain the best fit of the delays. In this latter case, the sinusoidal 
function has a period of 9 300 d, about half of that obtained using 
only the sinusoidal function. Moreover, the orbital period deriva¬ 
tive P = 1.44(7) x_10 -11 s/s ( compatible with P = 1.5(3) x 10~ u 
s/s obtained bv lHu et aPl2008l) is extremely high to be explained 
by a conservative mass transfer and loss of angular momentum 
from the binary system for gravitational radiation (see next sec¬ 
tion). This awkward result can be bypassed if the quadratic term 
is merely an approximation of a further sinusoidal function with 
a larger orbital period with respect to 9 300 d. 
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Folded period: 3000.6512 s 
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Fig. 6. Left and right plots show the folded RXTE/ASM light curve using LQS ephemeris (eq. [6jl and LSe ephemeris (eq. [8]) with P mo ,i = 18 600 
d, respectively. Each phase-bin is about 50 s. 


Under this assumption, the best fit obtained using the LQS 
ephemeris could be explained using a different scenario, where 
the quadratic term mimics the fundamental harmonic of a se¬ 
ries expansion whilst the sinusoidal term is the first harmonic. 
This seems also suggested by the best fit obtained using the LS 
function (eq.[3k since we obtain a modulation period, which is 
almost twice that obtained using the LQS function (eq.0. 

If we assume thatXB 1916-053 is part of a hierarchical triple 
system then the measured delays are also affected by the influ¬ 
ence of a third body. If the orbits of the third body and of the X- 
ray binary system around the common centre of mass are slightly 
elliptical then the delay A os (0 associated with the Doppler shift 
can be expressed as 


I e 

A-Ds(t) = A|sm(m f + w) + — [srn(2m f + m) - 3 sin(nj)] + 
e 2 

+ — [2 sin(3m, + m) - sin {m t + nr) cos(2 m, + 1)+ 

- 2 sin(m,) cos(ttr)] J, (7) 


where 

2n 

m, = -( t - t#) 

* mod 

is the mean anomaly; e is the eccentricity of the orbit; P m „d is the 
orbital period of both the X-ray binary system and the third body 
around the common centre of mass; m denotes the periastron an¬ 
gle; f (4 is the passage time at the periastron; and A — a sin i/c is 
the projected semi-major axis of the orbit, described by the cen¬ 
tre of mass of the X-ray binary system around the centre of mass 
of the triple system. We neglect third and higher order terms in 
Eq. |7] Limiting Eq.[7]to the first-order terms, it becomes the ex¬ 
pression shown bv lvan der Klis & Bonnet-Bidaudl(ll984l) . Then, 
we fitted the delays using 

y(t) = a + bt + A DS (t). 

Because the 27 available points do not cover a whole period, we 
arbitrarily fixed the value of P mod at 18 600, 17 100, and 20 100 


d, which are the best, lower, and upper value of the period ob¬ 
tained from the LQS ephemeris multiplied by a factor of two. 
The best-fit parameters are shown in Tab. [4] (columns 6, 7, and 
8 ). The ;p 2 (d.o.f.) are similar for the three adopted periods and 
the F-test probability of chance improvement with respect to LS 
function is 4.1 x 10~ 2 , 1.7 x 10~ 2 , and 0.9 x 10 2 for a P mo d value 
of 17 100, 18 600, and 20 100 d, respectively. In the following, 
we discuss the case of P mo d = 18 600 d. The best-fit function 
is indicated with a green curve in Fig. [3] Cleft panel). The corre¬ 
sponding residuals are shown in Fig. [3] (right panel, lower plot). 
The corresponding ephemeris (hereafter LSe ephemeris) is 

3 000.6512(6) 

T dip (N) = MJD(TDB) 50 123.010(3) +- V + 

86 4UU 

+A ds (N). (8) 

To verify the robustness of our results, we produced the 
folded light curves in the 3-5 and 5-12.2 keV energy bands of 
XB 1916-053 obtained from the All Sky Monitor (ASM) on 
board RXTE using the ephemerides shown above. We inferred 
those ephemerides using only pointing observations from which 
we obtained 27 points spanning from 1978 to 2014, whilst the 
RXTE/ASM light curves cover from 1996 Sep 01 to 2011 Oct 
31. We applied the barycentre corrections to the RXTE/ASM 
events. As a first step, we folded the RXTE/ASM light curves 
of XB 1916-053 using the LQ ephemeris reported by iHu et ail 
(l2008h and by us (Eq. QK adopting 60 phase-bins per period cor¬ 
responding to ~50 s per bin. The folded light curves and the 
corresponding hardness ratios (HRs) are shown in Fig. |4] None 
of the HR show an evident increase at phase zero as we would 
expect if the ephemerides well define the dip arrival times. This 
implies that those ephemerides do not correctly predict the dip 
arrival times contained in the RXTE/ASM light curve. Adopting 
the LQC ephemeris (eq. EK the maximum value of HR (that is 
2.8) is reached at phase 0.1 (see Fig. [5] left panels). Also in this 
case, the LQC ephemeris does not predict the dip arrival times in 
the ASM light curves of XB 1916-053. Using the LS ephemeris 
(Eq. [4j to fold the light curves, we obtained that the maximum 
value of HR is reached at phase zero and is close to 3.4 (see 
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Fig. 7. Folded RXTE/ASM light curve of XB 1916-053 selecting the events from SSC1 and SSC2. No energy filter is applied. Each phase-bin 
correspond to 75 s. Left panel: folded light curve using the LQS ephemeris (eq. 0. Right panel: folded light curve using the LSe ephemeris (eq. 

0 


Fig. 0 right panels). In contrast, with the LQS ephemeris (Eq. 
[6]» the maximum value of the HR falls in two phase-bins close 
to phase zero (see Fig. [6] left panels) and the maximum value of 
HR is 3.2, which is smaller than the value obtained with the LS 
ephemeris. Finally, we folded the RXTE/ASM light curves using 
the LSe ephemeris (eq. [8J. We show the folded light curves and 
the corresponding HR in Fig. [6] (right panel). In this last case the 
maximum value of the HR falls in only one phase bin at phase 
zero and the maximum value of the HR is about 4.5. 

We also folded the RXTE/ASM light curve (not filtered in 
energy) using the LQS and LSe ephemerides once we selected 
the events from the Scanning Shadow Cameras (SSCs) 1 and 
2. Adopting 40 phase-bins per period (that is each bin is 75 s), 
the folded light curves are very similar (see Fig. [7}, indicating 
that the two ephemerides are statistically equivalent. The dip is 
clearly observed at phase zero, the ASM count rate is reduced 
during the dip of 60% with respect to the persistent emission. Fi¬ 
nally, the goodness of the two ephemerides allows us to observe 
the presence of a secondary dip at phase 0.55, which is typically 
observed in several dipping sources (see iGrindlavi 19891 for XB 
1916-053). 


4. Discussion 

From the study of the 27 dip arrival times obtained from the 
pointed observations of XB 1916-053 and of the RXTE/ASM 
light curves, we find that the quadratic and cubic ephemerides 
do not correctly predict the dip arrival times on a long time span; 
whilst to well fit the delays, we need to use a function that con¬ 
tains at least linear and sinusoidal terms (LS ephemeris, see Eq. 
|4}. The addition of a quadratic term to the LS ephemeris (Eq. [6J 
gives a probability of chance improvement obtained with a F-test 
of 1.7 x 10 5 with respect to the LS ephemeris. Finally, using the 
ephemeris shown in Eq. [8] the probability of chance improve¬ 
ment, also with respect to the LS ephemeris, is 1.7 x 10 2 . The 
LQS and LSe ephemerides paint two different physical scenarios 
for XB 1916-053. In the first case the orbital period derivative of 


the X-ray binary system is P — 1.44(7) x 10~ u s/s and the ob¬ 
served delays associated with the dip arrival times are affected 
by a relatively low-amplitude (~ 130 s) sinusoidal modulation 
with a period close to 26 yr. In the second case the orbital pe¬ 
riod derivative is fixed to zero and the modulation of the delays 
is solely sinusoidal with an amplitude of ~ 550 s and an orbital 
period close to 51 yr. We explain in the following the sinusoidal 
modulation for both the scenarios, assuming the presence of a 
third body forming a hierarchical triple system with XB 1916- 
053, which alters the observed dip arrival times. 

We start by discussing the plausible values of the companion 
star mass Mt. We know that the companion star is a degenerate 
star and its radius /\S has to be equal to its Roche lobe radius R/ i 
since the binary system is in the Roche lobe overflow (RLOF) 
regime . Rearranging the Eq. 3.3.15 in IShapiro & Teukolskvl 
( 1983 1. the mass-radius relation for a degenerate star can be writ¬ 
ten as 


Ro IZ\ 5/i (M 2 

-^=0.04- —- 


R, 


\A) \M g 


-1/3 


= 0.0126 (1 +X) 5/3 m: 1/3 , 


where Z and A are the atomic number and the atomic weight 
of the matter composing the star, and where we assumed that 
the matter is only composed of hydrogen and helium. The factor 
Z/A is the average of Z/A for matter composed of hydrogen and 
helium, X is the fraction of hydrogen in the star and, finally, m 2 
is the companion star mass in units of solar mass. This equation 
has to be corrected for the thermal bloating factor /, which is the 
ratio of the companion star radius to the radius of a star with the 
same mass and composition, that is completely degenerate and 
supported only by the Fermi pressure of the electrons; then the 
factor / is > 1. The Roche lobe radius of the companion star can 
be written as 

R l1 = 0.46224 a (——— 

\ni\ +1112 

where a is the orbital separation of the binary system and m\ is 
the neutron star (NS) mass in unit of solar mass. We can write a 
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in terms of the orbital period P, m\, and m 2 , using Kepler’s third 
law. Combining the last two equations and Kepler’s third law, we 
obtain 

m 2 =0.0151 (1 +X) 5l2 f 12 . (9) 

iNelemans et al. 120061) , analysing the optical spectrum with the 
European Southern Observatory Very Large Telescope, detected 
a He-dominated accretion disc spectrum and suggested direct 
evidence for a helium donor. The authors found a good match 
with an LTE model consisting of pure helium plus overabundant 
nitrogen. For this reason, we assume X = 0 in the rest of the 
discussion. 

The bolometric X-ray flux of XB 1916-053 was estimated by 
several authors. iGallowav et alJ (120081) . analysing a RXTE/PCA 
observation of XB 1916-053, determined a persistent flux in the 
2.5-25 keV of (3.82 ± 0.04) x 10 “ 10 erg s ’ 1 cm- 2 . The authors 
corrected the flux for a bolometric factor <-■/,„/ = 1.37 ± 0.09 to 
estimate the bolometric flux in the 0.1-200 keV energy range, 
obtaining a_bolometric flux of (5.2 ± 0.3) x 10“ 1() erg s _1 cm“ 2 . 
Recently, IZhang et~akl (20140, analysing a Suzaku observation 
of XB 1916-053, found a value of Ft, 0 i in the 0.1-200 keV en¬ 
ergy range between 5.5 x 10 ~ 10 and 6.1 x 10 ~ 10 erg s“ ] cm -2 . 
Finally, analysing the persistent e miss ion of the source during a 
BeppoSAX observation, IChurch et al. (1998) estimated a value 
of Fboi in the 0.5-200 keV energy range of 6.2 x 10 “ 10 erg s -1 
cm -2 . Since the RXTE/ASM light curve of XB 1916-053 shows 
that the count rate of the source is almost constant over more 
than ten years, we adopt a conservative value for the bolometric 
flux of (5.5 ± 0.5) x 10 “ 10 erg s“' cm -2 . 

Th e distance d to the source was estimated by Gallowav et alJ 
(i2008h measuring the peak flux during the photospheric ra¬ 
dius^ expansion (PRE) in type-I X-ray bursts. Equation 8 in 
IGallowav et al. ( 2008 ) can be rewritten 


d = 8.32 


* pk,PRE 


3 X 10 ~ 8 


erg s“ 


-1/2 


1/2 


1 - 0.296 


m 1 


1/4 


fpRE) 

(1 +X)~ 1 / 2 kpc, 

( 10 ) 


where rp RE is the photospheric radius of the neutron star in units 
of 10 km and F p ^pr E is the flux at the peak of the type-I X-ray 
burst during the PRE. The authors measured F p kp RE = (2.9 ± 
0.4) x 1 0 “ 8 erg s“' crrT 2 and r PRE - 1.1 for XB 1916-053 and 
concluded that the distance to the source is d - 8.9 +1.3 kpc 
(adopting X — 0) for a NS mass of 1.4 M 0 . The X-ray luminosity 
can be expressed as L x = And 1 F/„,i, where we roughly assume 
that the emitted flux is isotropic. We obtain L x - 5.2 x 10 36 erg 
s _1 for a NS mass of 1.4 M e , whilst we find L x - 6.6 x 10 36 erg 
s _1 for a massive NS of 2.2 M Q . 

IRappaport et al.l d 19871 ) predicted the X-ray luminosity for 
highly compact binary systems under the reasonable hypothesis 
that the main mechanism to lose angular momentum is gravita¬ 
tional radiation. Combining the Eqs. 8 and 13 in their work, we 
obtain 

L, * , ^V W 4 ' 3 (1 + *)V // 3 erg S ± ( 11 ) 

where P m is the orbital period in units of minutes, [5 is the 
fraction of matter yielded by the companion star and accreted 
onto the NS, // is the efficiency for converting gravitational po¬ 
tential energy into X-ray emission, and a is the specific an¬ 
gular momentum carried away by the mass lost from the sys¬ 
tem, in units of 2na 2 /P or b, where a is the orbital separation (see 


IRappaport et al.lfl982 ). In Eq.[TT]we assume that the NS radius 
is 10 km. Using the orbital period value of 3 000.65 s, assuming 
r) — 1 and a conservative mass transfer scenario (fj = 1 ), we find 
that Lx — 1.1 X 10 35 / 3 erg s ' 1 and Lx — 2.3 x 10 35 / 3 erg s“ ] for 
a NS mass of 1.4 M e and 2.2 M e , respectively. Comparing the 
observed luminosity and the predicted luminosity, we estimate 
that / = 3.6 ±0.4 and / = 3.0+ 0.3 for a NS mass of 1.4 M Q and 

2.2 M 0 , respectively. Substituting the obtained values of / in eq. 
[9] we obtain that the companion star mass is M 2 = 0.10 ± 0.02 
M q and M 2 = 0.078 ±0.012 M G for a NS mass of 1.4 M Q and 2.2 
M 0 , respectively. The mass ratio q = M 2 /Mi of XB 1916-053 is 
between 0.036 ± 0.009 and 0.071 ± 0.009. 

iHu et al.l(i2008l) inferred the mass ratio of XB 1916-053 from 
the negative super-hump period and found q - 0.045, which is 
compatible with our estimated range of values of q. IChou et al] 
( 2001 1 estimated a value of q - 0.022 using the period of the 
apsidal precession of the accretion disc_of P prec — 3.9087(8) 
d. The value of q obtained bv IChou et al . ( 1200ill is outside the 
range that we find. 

To estimate the orbita l perio d derivative we use the eq. 11 
shown in IRappaport et al.1 (1 9871 ) that we rewrite as 

1 -l 4 5a(l Q -^) wfj ^ +I)5/2/3/2sS ~' (12) 

Using the value of P ~ 1.44 x 10 -11 s s _I (LQS ephemeris) 
and the orbital period value of 3 000.65 s, we find that the ther¬ 
mal bloating factor / is 40 and 32 for a NS mass of 1.4 and 2.2 
M 0 . These values of / are not physically plausible and suggest 
that, in a conservative mass transfer scenario, the value of the 
orbital period derivative cannot be that obtained from the LQS 
ephemeris. 

On the other hand, adopting an orbital period of 3 000.65 s 
and a factor / of 3.6 and 3.0 for a NS mass of 1.4 and 2.2 M G we 
find P = (3.9 ± 0.2) x 10 “ 13 s s“' and P = (3.98 ± 0.15) x 10 “ 13 
s s _I for a NS mass of 1.4 M G and 2.2 M 0 , respectively. The or¬ 
bital period derivative normalised to the orbital period is P/P - 

4.2 x 10 “ 9 yr 1 and weakly depends on the NS mass. We con¬ 
clude that the conservative mass transfer scenario with a thermal 
bloating factor of the companion star between three and four al¬ 
lows us to explain the discrepancy between the predicted and 
observed X-ray luminosity, but it does not solve the discrepancy 
between the predicted and measured orbital period derivative ob¬ 
tained from the LQS ephemeris. For this reason, we investigate 
the non-conservative mass transfer scenario. 

Combining the eqs. QT| and [T 2 ] we obtain 

y =* 3.38 x 10 51 /hiT7 „ 2 pf ,2 q erg s _1 . (13) 

Adopting L x ^ 5.2 x 10 36 erg s _1 , P = 1.44 x 10“ n s s _I , P = 
3 000.65 s and fixing q = 1, we find that y 8/ 3 ^ 2 = 0.191 for a 
NS mass of 1.4 M 0 . Since / > 1, we expect that more than 81% 
of the mass yielded by the companion star leaves the system. 
Furthermore, since the measured values of L x and P are positive, 
the term 1 - 1 ,5ar( 1 -/) in eqs. |TT] and [T 2 ] should be positive. 
Solving for a while taking /? < 0.191, we obtain that a < 0.823. 
Because a is in unit of 2na 2 /P or b, we find that the matter should 
leave the binary system from a distance d from the neutron star 
of d < a^ 2 a\ the point of ejection in unit of orbital separation is 
x = d/a < a 1 / 2 . In the rest of the discussion, we assume that the 
matter is ejected at the inner Lagrangian point x E \ of the binary 
system. We rewrite the eq.|TT]as function of / using the condition 
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Fig. 8. X-ray luminosity of XB 1916-053 in units of 10 36 erg s 1 versus 
the thermal bloating factor / of the companion star. The four curves 
correspond to different values of the NS mass: purple, green, light blue, 
and gold colours correspond to a NS mass of 1.4, 2, 2.1, and 2.2 M 0 , 
respectively. The peaks in the curves are at / ^ 1.5. 

/3f yl = 0.191. We find 



Fig. 9. Orbital period derivative of XB 1916-053 in units of 10 11 s s 1 
versus the thermal bloating factor / of the companion star. The brown 
and purple curves are obtained using a NS mass of 2.1 and 2.2 M 0 . 
The red and green lines indicate the best-fit value and the values at 68% 
confidence level of the orbital period derivative obtained from the LQS 
ephemeris. The purple curve is compatible at 1 cr with the measured 
orbital period derivative for / ^ 1.5. 


5.2 X 10 


,42 


1-1.5 X~Ll 2 (1-0.191 /- 3 /2) 


T| /3 P; l4/3 0.191/ 3/2 erg s -1 , 
(14) 


where xh is the position of the inner Lagrangian point in units 
of orbital separation. Using eq.[9]and a NS mass of 1.4 M 0 , xh 
can be written as a cubic function of / for values of the thermal 
bloating factor between 1 and 10. We find 

xh = 0.915 - 6.87 x 10~ 2 / + 6.61 x 10~ 3 / 2 - 2.88 x 10' 4 / 3 , 

with an accuracy of 2 x 10 3 . Combining the last equation and 
eq. [14j we infer the luminosity as function of /. We show L x in 
unit of 10 36 erg s _1 versus / for a NS mass of 1.4 M e (purple 
colour) in Fig. [8] Since the observed luminosity for a NS mass 
of 1.4 M 0 is larger than the predicted one for each value of f, 
also taking the corresponding error into account, we conclude 
that this specific non-conservative mass transfer scenario fails 
for a NS mass of 1.4 M Q . 

We repeat the same procedure for NS masses of 2, 2.1 and 
2.2 M q , finding that the predicted and observed luminosities are 
only compatible in the case in which the NS mass is > 2.2 M 0 . 
In this case ,we find that /?/ 3 / 2 = 0.154, a < 0.784 and 

xh = 0.927 - 6.02 x 10~ 2 / + 5.66 x 10~ 3 / 2 - 2.88 x 10" 4 / 3 

, with an accuracy of 2 x 10 5 . The luminosity for a NS mass of 
2.2 M 0 (gold colour) is shown in Fig. [8] Furthermore, we plot the 
orbital period derivative as function of / for a NS mass of 2.1 M Q 
(brown colour) and 2.2 M Q (purple colour) in Fig. [9] We note that 
only for a NS mass of 2.2 M e the predicted and measured P are 
compatible for / ^ 1.5. We conclude that this non-conservative 
mass transfer scenario predicts the observed values of luminosity 
and orbital period derivative only for NS masses larger than 2.2 
M 0 . For a NS mass of 2.2 M Q , the companion star has a mass of 
0.028 M 0 and ft is close to 0.084, which is more than 90% of the 
matter, yielded from the companion star, that leaves the binary 
system. 

In this scenario, we suggest that XB 1916-053 could 
be considered as a possible progenitor of the ultra-compact 
"Black Widow" pulsars with very low-mass companions. 


[Benvenuto et al.l (1201 2|) proposed that a binary system with an 
initial orbital period of 0.8 d, composed of a 1.4 M Q NS and 
a companion star mass of 2 M 0 , evolves in ~6.5 Gyr forming 
a binary system that well fits the known orbital parameters of 
the black widow millisecond pulsar PSR J1719-1438. We note 
that the same evolutive path fits the orbital parameters of XB 
1916-053 at ~ 5 Gyr from the initial time. At 5 Gyr, the pre¬ 
dicted orbital period is 0.035 d, the predicted companion star 
mass is 0.03 M e , the NS mass is slightly larger than 2.2 M G 
(Benvenuto, private communication) and the companion star is 
helium dominated. These values are very similar to those of 
XB 1916-053 shown in this work for a non-conservative mass 
transfer scenario, although a discrepancy between our estima¬ 
tion of Mi ~ 4.1 x 10 ~ 9 M 0 yr and the value suggested by 
iBenvenuto et al.l (12012i) at 5 Gyr (~ 10 ~ 10 M G yr *) is present. 
Furthermore, we note that as the spin period of PSR J1719-1438 
is 5.7 ms (see lBailes et al.ll201 ll and references therein) the spin 
period of the NS in XB 1916-053 could also be extremely short. 
Indeed, iGallowav et al.l (1200lb interpreted the asymptotic fre¬ 
quency of the coherent burst oscillations in terms of a decoupled 
surface burning layer and suggested that the NS could have a 
spin period around 3.7 ms. 

Nevertheless, we note that our solution for a non¬ 
conservative mass transfer scenario is not supported by a ro¬ 
bust physical mechanism to explain the large quantity of mat¬ 
ter ejected from the inner Lagrangian point. To date, only two 
physical mechanisms are known to be able to eject the trans¬ 
ferred matter partially (or totally). The first mechanism predicts 
that when a super-Eddington mass transfer occurs, the X-ray lu¬ 
minosity has to be at the Eddington limit. Then, the radiation 
pressure from the compact object pushes away part of the trans¬ 
ferred matter from the binary system. This mechanism was re¬ 
cently invoked to explain the large orbital period derivative mea- 
sured in the accretio n disc corona (ADC ) so urce XI822 - 371 b y 
iBurderi et al.l (I2010I) . Ilaria et aD (1201 3b . and llaria et akl (1201 5b . 
However, this mechanism cannot be applied in the case of XB 
1916-053 because type-I X-ray bursts are observed in the light 
curve of the source (see e.g. Fig. 0 , whilst the stable burning sets 
in at high accretion rate values that are comparable to the Ed¬ 
dington limit (see Bildstenll2000l and references therein). Con- 
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Table 5. Best-fit parameters of the delays assuming the presence of the 
third body in eccentric orbit and taking a quadratic term c = 5 x 10 -7 
s/d 2 into account. 


Parameters 

Pmod= 17 100 d 

P mod =18 600d 

P mod = 20100 d 

a (s) 

180 + 332 

21 ±307 

-27 ± 285 

Z?(xl 0~ 3 s/d) 

2 + 20 

2 + 19 

4 + 21 

A(s) 

506 ± 46 

534 ± 43 

562 ± 43 

e 

0.26 ± 0.20 

0.28 ±0.15 

0.32 + 0.13 

m (deg) 

198 + 27 

213 + 28 

219 + 27 

U (d) 

-3 594+ 1 129 

-3 036+ 1 131 

-2 825 ± 1 042 

T 2 (d.o.f.) 

51.3(21) 

47.9(21) 

45.5(21) 

F-test prob. 

3.5 x IQ ” 2 

1.5 x 10 ' 2 

0.8 x 10~ 2 


Note — The reported errors are at 68% confidence level. The F-test 
probability is estimated with respect to the^ 2 value of the LS ephemeris 
(the fourth column of Tab. 0. 


sequently, the mass transfer rate cannot be super-Eddington and 
this mechanism cannot justify a non-conservative mass transfer 
scenario. The second mechanism supposes that the X-ray binary 
system is a transient source and during the X-ray quiescence 
it is ejecting the transferred matter from the inner Lagrangian 
point due to the radiation pressure of the magneto-dipole rota¬ 
tor emissioi^Thisjnechanism, which we call radio eject i on af- 
ter lBurderi et al.l (12001 1) , was proposed bv lDi Salvo et al.1 (f2008) 
to explain the large orbital period derivative measured in SAX 
J1808.4-3658. However, this mechanism also fails to explain our 
results because XB 1916-053 is a persistent X-ray source. 

Finally, we discuss the sinusoidal modulation observed in the 
LQS and LSe ephemerides. If we assume a conservative mass 
transfer scenario, the predicted orbital period derivative is close 
to 4 x 10 - 13 s s 1 independent of the NS mass. Then we added a 
quadratic term to the LSe ephemeris to take the predicted value 
into account. We fitted again the delays using the relation 

y(t) - a + bt + ct 2 + A ds ( t ), 

where the term c is fixed to 5 x 10 7 s/d 2 . The fit parameters 
are reported in Tab. [5] We note that the addition of the quadratic 
term does not significantly change the best-fit parameters. 

An explanation of the sinusoidal modulation obtained from 
the LSe ephemeris could be the presence of a third body gravi¬ 
tationally bound to the X-ray binary system. Assuming the exis¬ 
tence of a third body of mass M 3 , the binary system XB 1916- 
053 orbits around the new centre of mass (CM) of the triple sys¬ 
tem. The distance of XB 1916-053 from the new CM is given 
by a x = cibi,, sin i = A c, where i is the inclination angle of the 
orbit with respect to the line of sight, A is the amplitude of the 
sinusoidal function obtained from the ephemeris of eq. [ 8 ] and c 
is the light speed. We obtained a x = (1.60 ± 0.13) x 10 13 cm for 
Pmod = 18 600 d. We can write the mass function of the triple 
system as 

M 3 sin/ / 47 t 2 \^ 3 a x 

(m 3 + M bin y-p = j 

mod 

where M 3 is the third body mass, M&,-„ the binary system mass, 
and finally, P mo d is the orbital period of XB 1916-053 around the 
CM of the triple system. Substituting the values of M;„„, P mo d, 
a x , and assuming an inclination angle for the source of 70°, we 
find that m 3 is ~ 0.10 M e and ~ 0.14 M e for a NS mass of 1.4 
M 0 and 2.2 M 0 , respectively. We used also P mo d of 17 100 and 
20 100 d finding that the values of m 3 are substantially indepen¬ 
dent of the value of P mo d- 


For a non-conservative mass transfer scenario, we discuss 
the sinusoidal modulation obtained from the LQS ephemeris as¬ 
suming a NS mass of 2.2 M 0 . In this case we find that a x = 
(3.9 + 0.5) x 10 12 cm and m 3 ~ 0.055 M Q for an inclination angle 
of 70°. 


5. Conclusions 

We have systematically analysed all the historically reported X- 
ray light curves of XB 1916-053, which span 37 years. We find 
that the previously suggested quadratic ephemeris for this source 
no longer fits the dip arrival times. 

We studied the conservative mass transfer scenario of the 
system, finding that the thermal bloating factor of the degenerate 
companion star is 3.6 and 3 for a NS mass of 1.4 and 2.2 M Q . In 
this scenario, the predicted and observed luminosity are compat¬ 
ible (~5-7 xlO 36 erg s 1 ), although the orbital period derivative 
is a factor of 40 smaller than the value of 1.44 x 10~ n s s _1 ob¬ 
tained fitting the delays with a quadratic plus a sinusoidal func¬ 
tion (LQS ephemeris). If the conservative mass transfer scenario 
is correct, we conclude that the modulation of the delays asso¬ 
ciated with the dip arrivals time are solely due to a sinusoidal 
modulation caused by a third body orbiting around the binary 
system. In this case we estimate the third body mass is 0.10 and 
0.14 M 0 for NS masses of 1.4 and 2.2 M e , respectively. The or¬ 
bital period of the third body around XB 1916-053 is close to 55 
yr and the orbit shows an eccentricity e = 0.28 ± 0.15. 

In a non-conservative mass transfer scenario where the mass 
is ejected away from the inner Lagrangian point, we find that the 
observed luminosity and the orbital period derivative obtained 
from the LQS ephemeris are possible only from a NS mass > 2.2 
M q . In this case we obtain that the thermal bloating factor of the 
degenerate companion star is f - 1.5, the companion star mass 
is 0.028 M 0 , and the fraction of matter yielded by the compan¬ 
ion star and accreting onto the NS is f5 — 0.084. In this scenario, 
the sinusoidal modulation of the delays can be explained by the 
presence of a third body orbiting around XB 1916-053 with an 
period of 26 yr. We find that the third body mass is 0.055 M e . Fi¬ 
nally, if the non-conservative mass transfer scenario is valid, we 
suggest that XB 1916-053 and the ultra-compact black widow 
system PSR J1719-1438 could be two different stages of the 
same evolutive path discussed by [Benvenuto et alJ ( 2012 ’). If it 
is true, then the age of XB 1916-053 is close to 5 Gyr, whilst 
PSR J1719-1438 is ~6.5 Gyr old. 
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